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We study here the effect of a varying G on the evolution of the early Universe and, in particular, on 
primordial nucleosynthesis. This variation of G is modelled using the Brans-Dicke theory as well as 
a more general class of scalar-tensor theories. Modified nucleosynthesis codes are used to investigate 
this effect and the results obtained are used to constrain the parameters of the theories. We extend 
previous studies of primordial nucleosynthesis in scalar-tensor theories by including effects which can 
cause a slow variation of G during radiation domination and by including a late-time accelerating 
phase to the Universe's history. We include a brief discussion on the epoch of matter-radiation 
equality in Brans-Dicke theory, which is also of interest for determining the positions of the cosmic 
microwave background power-spectrum peaks. 
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I. INTRODUCTION 

The possibility of the physical "constants" taking dif- 
ferent values at different times in the Universe's history 
has recently received much attention with the apparent 
observation that the fine structure constant had a dif- 
ferent value in the distant past pj. These observations, 
if correctly interpreted , motivate further study of the 
effects of any variation on other physical constants. The 
scalar-tensor (ST) theories of gravity provide us with a 
self-consistent way of modelling a possible variation in 
Newton's constant, G. In these theories the metric ten- 
sor field of general relativity (GR) is no longer the only 
field mediating the gravitational interaction; there is an 
additional scalar component. Test particles in these the- 
ories still follow geodesies of the metric but the manner 
in which the metric is generated from the sources is al- 
tered by the presence of the scalar field. Such theories 
were first considered by Jordan in 1949 Q and later de- 
veloped by Brans and Dicke in 1961 for the case of 
a constant coupling parameter. The more general case 
of a dynamic coupling was considered by Bergmann , 
Nordtvedt @ and Wagoner [j| in the 1960s and 70s. 
Solution-generating procedures for these theories were 
found by Barrow and Mimoso Barrow and Parsons 
[HI and Barrow [12. 

Using these theories we investigate the earliest well 
understood physical process, primordial nucleosynthesis. 
Previous studies on this subject have been carried out 
by a number of authors. In particular, the Brans-Dicke 
(BD) theory has been especially well studied in this con- 
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text W, for example, Casas, Garcia-Bellido and Quiros 
[I3 , uM an d Serna, Dominguez-Tenreiro and Yepes ^| . 
The more general class of ST theories has also been well 
studied, most notably by Serna and Alimi |16| . Santi- 
ago, Kalligas and Wagoner 0] and Damour and Pichon 
[fj| . More recently, a number of studies have been per- 
formed investigating the effect of a nonminimally cou- 
pled quintessence field on primordial nucleosynthesis (see 
e.g. [l3)|2fj)- Although these studies are very detailed 
they all make the simplifying assumption of a constant G 
during the radiation-dominated phase of the Universe's 
history (with the exception of |15|. who numerically in- 
vestigate the effect of an early scalar-dominated phase 
on the BD theory, and 0] who use the idea of a "kick" 
on the scalar field during electron-positron annihilation). 
We relax this assumption and investigate the effects of 
entering the radiation-dominated phase with a noncon- 
stant G. The constraints we impose upon the variation 
of G during primordial nucleosynthesis are then used to 
constrain the parameters of the theory. In carrying out 
this study we consider the more general class of ST the- 
ories, paying particular attention to the BD theory. 



II. SCALAR-TENSOR GRAVITY 

The ST theories of gravity are described by the action 




+ WnL m ). (1) 

Here, is the scalar field, R = is the Ricci scalar 
and L m is the Lagrangian describing matter fields in the 
space-time. The action above can be extremized with 
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respect to g^ to give the field equations 



III. COSMOLOGICAL SOLUTIONS 
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and with respect to (f> to give scalar-field propagation 
equation 
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where a prime denotes differentiation with respect to (f>, 
T^ u is the energy-momentum tensor of the matter fields, 
defined in the normal way, and T = T M M . In the limit 
u) — > oo and uj' /uj 3 — > the ST theories reduce to GR. 

A useful property of these theories is that they can 
be conformally transformed to a frame where the field 
equations take the same form as in GR, with a nonmin- 
imally coupled, massless scalar field. We refer to these 
two frames as the Jordan (J) frame and the Einstein (E) 
frame, respectively. Under the conformal transformation 
9iiu = A 2 {<j))g^ Ul where A 2 {<j)) = 1/0, the Einstein - 
Hilbert action for GR containing a scalar field, ijj, and 
other matter fields, *ff m , 
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1 

16tt 



■S m [^ m ,A 2 g llv ]. (4) 



is equivalent to the ST action (I}, if we make the defini- 
tions r = \nA{<j)), \f^Ka = dT/dijj and the identification 
or 2 = 2cu((f>) + 3. Here, symbols with bars refer to quan- 
tities in the E frame and symbols without bars refer to 
quantities in the J frame. 

By extremizing the action with respect to g^ we 
get the E-frame field equations 

Rpv = -8tt(T^ - ^g^f + V'.mVv) 

and by extremizing with respect to ip we get the E frame 
propagation equation 

Dtp = -Vl^af, 

where we have defined the energy-momentum tensor in 
the E frame with respect to g^ so that f^ v = A 6 T^. 

It is worth noting that whilst the J-frame energy- 
momentum tensor is always covariantly conserved, 
T^ v . v = 0, its counterpart in the E frame is not, T'" y . ly = 

^/inaTtp,^ . By insisting that in a physical frame the 
energy-momentum of the matter fields should always be 
conserved we are forced to identify the J frame as the 
physical one. Any results derived in the E frame must be 
transformed appropriately, if they are to be considered 
as physical. 



A. Brans— Dicke theory, uj ^constant 

Consider the homogeneous and isotropic space-times 
described by the J-frame Friedmann-Robertson- Walker 
(FRW) line-element 



ds 2 = dt 2 
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where k = 0, ±1 is the curvature parameter. 

Substituting into J2J) and J2J gives, for a perfect 
fluid, the J-frame Friedmann equations 



a fa 
2-+ - 

a \a 



u> 4> 2 a 4> 



8tt 



+ 2 -T + T = --T P-—, (6) 

2 Z a <p <p <p a z 



8tt 



and 



d\ a (f> uj (f> 2 k 
a <p 6 <j) 2 a 2 ' 



= 8tt (p - 3p) _ g a0 
<\> 4> (2w + 3) o0' 



(7) 



(8) 



where overdots denote differentiation with respect to t. 
From T^ v - V = we also obtain the perfect fluid conser- 
vation equation 



p + 3-(p + p)=0. 
a 



1. Radiation domination 

We consider first the case of a flat radiation-dominated 
universe. Assuming the equation of state p — |p and 
defining the conformal time variable, 77, by adn = dt, the 
equations (jSJ, Q and (|SJ integrate to give 
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a{rj) rf + (2tj 2 + 3771)77 + r?f + 3771772 - \unf x 
<f/(rj) 37/i 
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(73(77) T] 2 + (2772 + 3771)77 + r}\ + 3?7i 772 - §^77? 

where 771 and 772 are integration constants. 

For uj > —3/2 the solutions to these equations are 
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i(r,)=a 1 (r l + r )+ ) V^r, + r?_) 2 V^F (n) 



4>kn)=<h{v + v+) 2 v^(v + v-) (12) 



where ij± = 772 + §771 ± §?7iy 1 + a i and <j>\ are in- 
tegration constants, and 87rp rO /30iai = 1 (subscript 
indicates a quantity measured at the present day and we 
rescale so that ao = 1, throughout). 
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FIG. 1: The evolution of a and cj> as functions of rj for u> > 
-§ with \rn \ = 10 -12 C?eV _1 , u = 40000, and r) 2 set so that 
a(0) = 0. The solid line corresponds to rji > and the dashed 
line to r\\ < 0. 



FIG. 2: The evolution of a and <j) as functions of rj for ui < 
-§ with \rji\ = 10~ 13 Gel/ _1 , u = 40000, and n 2 set so that 
a'(0) = 0. The solid line corresponds to r/i > and the dashed 
line to rji < 0. 



For uj < —3/2 wc find 



a(rf) — ai exp 



-1 -iV + V- 

x exp | — , tan 
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where 77+ — 2 , ' 1 V 3 
87r/9 r o/30ia^ = 1. 

In these solutions 771 determines the evolution of the 
scalar field during radiation domination, and is a physi- 
cally interesting quantity; 7/2 sets the origin of the con- 
formal time coordinate. The evolution of a and <b, for 



oj > —3/2 and to < —3/2 , is shown in figures ^ and [21 
respectively. A value of u> = 40000 is chosen, in agree- 
ment with recent solar system observations [22j. 

For oj > —3/2 (lo < —3/2) we see that the scale factor 
here undergoes an initial period of rapid (slow) expan- 
sion and at late times is attracted towards the solution 
(2(77) oc 77, or a(t) cx t?. Similarly, <f> can be seen to be 
changing rapidly at early times and slowly at late times. 
We attribute these two different behaviours, at early and 
at late times, to periods of free scalar-field domination 
and radiation domination, respectively. In fact, setting 
771 = in @ and (|10fl . we remove the scalar-dominated 
period and gain the power-law exact solutions 0(77) cx 77 
and 4> — constant. These are the "Machian" solutions 
usually considered in the literature. Allowing 771 to be 
nonzero we will have a nonconstant 0, and hence G, dur- 
ing primordial nucleosynthesis. If p r o — is chosen then 
these solutions become vacuum ones that are driven by 
the cf> field. For k — these solutions do not have GR 
counterparts. 

We see from (|14H . and figure [21 that for uj < —3/2 the 
initial singularity is avoided; for a more detailed discus- 
sion of this effect see 1231 . 
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2. Matter domination 



and 



When considering a matter-dominated universe we 
could proceed as above and determine a set of general 
solutions to JfJJl, and © that at early times are de- 
scribed by free scalar-field domination and at late times 
by matter domination (see |2lj|). but for our purposes 
this is unnecessary. For a realistic universe we require a 
period of radiation domination during which primordial 
nucleosynthesis can occur. If the scalar-field-dominated 
period of the Universe's history were to impinge upon the 
usual matter-dominated period then we would effectively 
lose the radiation-dominated era. For this reason it is 
sufficient to ignore the free scalar component of the gen- 
eral solution and consider only the Machian component. 
This is equivalent to imposing the condition tfia 3 — > as 
a — ► 0. With this additional constraint the solutions to 
©, 10 and @, for k = and p = 0, are given by 



if) + 3—tjj = —V4na(p — 3p) 



(19) 



a(t) = a*i 4 + 3 " and (f>(t) = <t>*t 4 + 3 «> (15) 



where 



8*_PrnU = 2(3 + 2hQ 

30, al ~ 3(4 + 3u) 

and a* and are constants. These solutions can be seen 
to approach their GR counterparts as ui — > oo. 



3. Vacuum, domination 

Similarly, for a vacuum (p = —p) dominated period of 
expansion we can impose the condition 4>a 3 — > as a — > 
to get the power law exact solutions 24] 



a(t) = a^ +UJ and tf>(t) = t r 



(16) 



where 



Stt 1 

^Pao = ^(3 + 2^)(5 + 6w). 

and a-f and </>j are constants. Note that this vacuum 
stress does not produce a de Sitter metric as in GR. 



B. Dynamically— coupled theories, u> = co((j>) 

In order to evaluate more general scalar-tensor the- 
ories it is convenient to work in the E frame. In this 
frame, substituting the FRW line-element into the field 
equations gives 



1 
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where over-dots here denote differentiation with respect 
to t and a(i) = A(cp)a(t) is the scale-factor in the E 
frame. Defining N — ln(a/ao) Damour and Nordtvedt 
H| write |T7 |) . and <P|l as 

w^/' +(2 -^h*'=- (i - 3 ->^ (20) 



where e = 3fc/87rpa 2 , p = (7— l)p and ' denotes differen- 
tiation with respect to N . 

Now we consider a coupling parameter of the form 

ro/0 = r(Voo) + oo(v(t) - Voo) + \f3(m - Voo) 2 

so that 



a(t) = a + —={ij)(t) - ipoo). 

V47T 



(21) 



This is of the same form as chosen by San tiag o, Kalligas 
and Wagoner 01 an( i Damour and Pichon |18| . Santiago, 
Kalligas and Wagoner arrive at this form of a by assum- 
ing the evolution of the Universe has been close to the GR 
solutions throughout the period from primordial nucle- 
osynthesis to the present. They therefore consider them- 
selves justified in performing a Taylor expansion about 
the asymptotic value of ip and discarding terms of sec- 
ond order or higher. Damour and Pichon consider T to 
be a potential down which ip runs. They assume that 
a particle near a minimum of a potential experiences a 
generically parabolic form for that potential. This leads 
them to consider a quadratic form for T which gives, on 
differentiation, a as above. These two lines of reasoning 
are, of course, equivalent as the parabolic form of a par- 
ticle near its minimum of potential can be found using 
a Taylor series. Theories with this form of a belong to 
the attractor class which approach GR at late times, if 
we impose the additional condition «o = 0. 



1. Radiation domination 

For the case of a radiation-dominated flat universe, the 
general solutions of l(T7|) . (fTHfl and (|T$|) with the choice 
PTJl. are, for ui > -3/2: 

a 2 (v) = ~~7r~( r i - m){v -m + M) 



V-V2 
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°>{v) = — — m\ +^r~ (v-m) 
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Here, 771 determines the evolution of the scalar held, tpi 
is the value it approaches asymptotically and 772 sets 
the origin of the conformal time coordinate defined by 
a(r])drj = di. For u> > —3/2, 771 6 1Z such that ip € TZ 
whilst for uj < —3/2, 771 E X such that ip E X and weK. 
If the idea of an imaginary scalar held seems strange 
the reader is reminded that the E frame is considered 
to be unphysical: all observable quantities measured in 
the physical J frame remain real at all times. 
These solutions in the J frame are 



4>(ri) = exp (-f3(ip(rj) - V>oo) 2 ) 



and 



a 2 (r?) 



(22) 



(23) 



They exhibit the same features as their BD counterparts: 
at early times there is a period of free-scalar-held dom- 
ination, and at late times they approach oat' and 
6 =constant. 



2. Matter domination 

Solutions during the matter-dominated era are diffi- 
cult to hnd because the energy-momentum tensor for 
the matter held is not conserved. However, it is possible 
to obtain an evolution equation for ip. Using 12U|) . we get 



-^M^" + (l-~^' = -^-VO. (24) 



(3-47Ti// 2 ) 



4tt 



For a flat universe e = and we can solve ()24|l by mak- 
ing the simplifying assumption tp' 2 <C 3. This gives the 
solution 



+ Be-^-V^W)* (25) 
where A and B are constants of integration. 



3. Vacuum, domination 

Similarly, for the case of a vacuum-dominated universe 
the evolution of ip can be approximated using (|20|1 to 
obtain 



+ De -Hi-V^W)N (26) 
where C and D are constants of integration. 



IV. PRIMORDIAL NUCLEOSYNTHESIS 

A. Modelling the form of G(t) 

If primordial nucleosynthesis were to occur during the 
scalar-dominated period then the very different expan- 
sion rate would have disastrous consequences for the 
light-element abundances (see e.g. Therefore we 

limit our study to times at which the scale factor can 
be approximated by a form that is close to a(rj) oc 77, 
and so primordial nucleosynthesis can safely be described 
as occurring during radiation domination. Performing a 
power-series expansion of the solutions l|llf> . (|12|l . i|13|) 
and l|14fl in 771/(77 + 772) we hnd 

a{j]) = 01(77 + t?2 + 3?7i) +0(rjf) (27) 

3r h "\ , rM 2s 



0(77) 



1 



77 + 772 



(28) 



for both uj > —3/2 and lu < —3/2. We can then set 
the origin of the 77 coordinate such that a(0) — with 
the choice 772 = — 3?7i. The solutions (|2T| and f2"%)l then 
become, in terms of the time coordinate t, 

a(t) =ait* +0(?7?) 

*(i)=* (l + ^)+°<rf). 

where 02 = — 3?7iai and the origin of t has been chosen 
to coincide with the origin of 77. 

Similarly, expanding the solutions 12211 and l|23() in 
771/(77 + 772) we hnd that, for both uj > —3/2 and 
u) < -3/2, 

a(t) = ai*i +0{rjl) 
a(t) 

where 



</3(t) = 4> x 1 + + 0{rf x ) 



ax = v/87rp r o/3exp(/3(7/)i - t/> oc ) 2 /2) 

02 = V 2p r0 (3r)i (ipi - t/>oo) exp(/3(7/'i - ip a 

</>i = exp(-/3(V>i - tAoo) 2 ) 



2 /2) 



and 772 has been set so that a(0) = 0. 

In the limit <f> — > we can see from equations Q 
that the standard GR Friedmann equations are recov- 
ered, with a different value of Newton's constant given 
by 



G(f) = W) 



Gi 



o(t) 



a(t) + a 2 



(29) 



where Gi — l/<f>i- We conclude that the solutions found 
above correspond to a situation that can be described 
using the GR Friedmann equations with a different, and 
adiabatically changing, value of G. This is just the situa- 
tion considered in recent work by Bambi, Giannotti and 
Villante |26| . but we now have an explicit form for the 
evolution of G(t) derived from ST gravity theory. 
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FIG. 3: The upper and lower bounds on Gi/Gn as a function of ai are shown by the solid lines. The shaded regions correspond 
to the allowed parameter space and the dashed lines show Gi = Gn- 



B. The effect of a non— constant G 



and the 4 He abundance estimated by Barger et al. [2E 



The time at which weak interactions freeze out in the 
early Universe is determined by equality between the rate 
of the relevant weak interactions and the Hubble rate. 
When the weak interaction rate is the greater then the 
ratio of neutrons to protons tracks its equilibrium value, 



Y P = 0.238 ± 0.005 to la, 



(31) 



n/p = exp(— (to,, 



,)/T) where m T . 



and rrip are the 
neutron and proton masses. If, however, the Hubble rate 
is greater than the weak-interaction rate then the ratio 
of neutrons to protons is effectively "frozen-in" , and /3 
decay is the only weak process that still operates with 
any efficiency. This will be the case until the onset of 
deuterium formation, at which time the neutrons become 
bound and (3 decay ceases. The onset of deuterium for- 
mation is primarily determined by the photon to baryon 
ratio, ry 7 , which inhibits the formation of deuterium nu- 
clei until the critical temperature for photodissociation is 
past. As the vast majority of neutrons end up in 4 He the 
primordial abundance of this element is influenced most 
significantly by the number of neutrons at the onset of 
deuterium formation, which is most sensitive to the tem- 
perature of weak-interaction freeze-out, and hence the 
Hubble rate, and so G, at this time. Conversely, the pri- 
mordial abundances of the other light elements are most 
sensitive to the temperature of deuterium formation, and 
hence ?y 7 , when nuclear reactions occur and the light el- 
ements form. (See |26| for a more detailed discussion of 
these points). 

Using the simple forms of a(t) and G(t) derived above 
we use a modified version of the Kawano code |27j | to in- 
vestigate the effect of this variation of G on primordial 
nucleosynthesis directly. We use the deuterium abun- 
dance estimated by Kirkman et al. |2^| 



log(f 



= -4.556 ±0.064 tola, 



(30) 



to create a parameter space in (r/j, G±, 02). The three- 
dimensional 95% x 2 confidence region is projected into 
the G\, a-i plane to give figure |3J where three species of 
light neutrinos and a neutron mean lifetime of r = 885.7 
seconds have been assumed. 

The apparent disfavoring of the value of G\ much 
greater than Gn in figure is due to the observational 
limits we have adopted for Yp (equation l|31|0 . which 
suggest a helium abundance that is already uncomfort- 
ably low compared to the theoretical prediction for Yp 
derived from the baryon density corresponding to equa- 
tion |201 We expect these constraints to be updated by 
future observations (see e.g. 30] ); this will require a cor- 
responding update of any work, such as this, that seeks 
to use these constraints to impose limits on physical pro- 
cesses occuring during primordial nucleosynthesis. 

These limits on the parameters G\ and ai can be used 
to construct plots showing the explicit evolution of G(t) 
for various limiting combinations of the two parameters, 
this is done in figure 21 It is interesting to note that 
in both of these plots the lines corresponding to differ- 
ent values of (12 all appear to cross at approximately the 
same point, log a ~ —9.4. This confirms our earlier dis- 
cussion, and the results of [2(j, that the 4 He abundance is 
mostly only sensitive to the value of G at the time when 
the weak interactions freeze out. In reality, this freeze- 
out happens over a finite time interval, but from figure 
01 we see that it is a good approximation to consider it 
happening instantaneously - where the lines cross. To 
a reasonable accuracy one could then take G through- 
out primordial nucleosynthesis to be its value during the 
freeze-out process. 
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FIG. 4: These plots show the explicit evolution of G for some limiting values of Gi, as taken from figure^ for various values 
of a 2 . Solid lines correspond to |a 2 j = 10 -11 , dashed lines correspond to |a 2 | = ^q-ii-155 an( ^ dotted lines correspond to 
|a 2 | = 10- 12 . 



V. CONSTRAINING THE THEORY 

Using the results in the previous section it is possible to 
constrain the underlying ST gravitational theory. This 
is done separately for the BD theory and for the more 
general ST theories. For each theory we consider the 
case of a universe containing matter and radiation only 
and then a universe containing matter, radiation and a 
non-zero vacuum energy, with equation of state p = —p. 



A. Brans— Dicke theory 

1. Universe containing matter and radiation 

Using equations (|15|l we can write the ratio of G at 
matter-radiation equality, G eg , to its present-day labo- 
ratory value, Gat, as 

(2w + 3)^ (2^ + 3) / a \ 



Gjv (2cj + 4) cj> eq (2w + 4) V a eq 



(2cj + 3) f 32 1 
(2lj + 4) U + ^) (32j 



where we have used the expression for G in the weak-field 
limit i3 



G(0 = 



(2cj + 4) 1 
(2w + 3)^' 



for Gat, whilst using the expression G(<f>) — l/4> for G eq , 
as in 129( 1 . In reality, the laboratory value of G today 
is not equal to that in the background cosmology |32|| . 
although we take it to be so here for simplicity. 

We now proceed by calculating 1 + z eq in the Brans- 
Dicke cosmology, following Liddle, Mazumdar and Bar- 
row [33|. As T^ v . v = in the J frame we have that 



p r oc a 4 and p m oc a a . For a universe containing mat- 
ter and radiation only this gives 



1 + Zeq 



a, 



eq 



PmQ Preq 
PrO Pmeq 



PmO 
PrO 

, when entropy 



as well as the usual relation T oc a~ 
increase is neglected. 

Photons and neutrinos both contribute to the value of 
p r o, the present-day energy-density of radiation. From 



T, 



70 



2.728 ± 0.00AK 



we get p 7 o 



4.66 x 



10 34 g cm 3 and using the well known result T v0 — 
(4/ll) 1 / 3 T 7 o, and assuming three families of light neu- 
trinos, we also have p u $ = 0.68p 7 o- This gives the 
total present-day radiation density as p r o — 7.84 x 
10 _34 g cm~ 3 . Now, recalling our assumption of spatial 
flatness and J7J|, we can write 



PtotO — PmO + PrO 



3# 2 (4 



-3w)(4 + 2w) 



BttGat 6(1 + uj) 2 



For G N = 6.673 x 1CT 11 Nm 2 kg- 2 , H = 
lOO/i kms' 1 Mpc~ x , and the value of p r o above, we have 



1 



= 2.39 x 10% 



(4 + 3w)(4+ 2oj) 



6(1 



1 



= 2.39 x 10 4 h 2 1 



(33) This correction to 1 + z eq has direct observational con- 



sequences in the power spectrum of cosmic microwave 
background (CMB) perturbations. After l + z eq the sub- 
horizon scale perturbations, that were previously effec- 
tively frozen, are allowed to grow. Changing the value 
of 1 + z eq therefore causes a shift in the power-spectrum 



eaks, which is potentially observable (see |3^, [3^|, [3 



peaks, wi 

[33 , [38| and [39] for a more detailed discussion of the ef- 
fect of a varying G on CMB formation) . For our purposes 
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this modified expression for 1 + z eq can then be substi- 
tuted into (|32[1 to give an equation in terms of u>, G eq 
and h. Assuming a value of h = 0.7 and that G eq = G\, 
as appears a very good approximation for the models 
above, we can use our bounds on G\ in terms of a 2 to 
create an allowed parameter space in the {01,0,2) plane. 
This is shown in figure 

We remind the reader that the a large u> corresponds to 
a slowly varying Machian component of 0(1 + z), as can 
be seen directly from equations l|15|l . i|16|) and the late- 
time limits of (JT1J, flit ) CSl an d 03} ■ The parameter a 2 
determines the evolution of the free component of 0, as 
can be seen from the early time limits of i|ll|) , l|12l) , 
and JHJ. In the limits u> — ► 00 and a 2 — ► the Machian 
and free components of 0(1 + z) both become constant, 
respectively, resulting in a constant G(l + z). The con- 
straints imposed upon ui and 0,2 in figure |3] therefore cor- 
respond to constraints upon the evolution of G(l + z) in 
this theory, valid both during priomordial nucleosynthe- 
sis and at other cosmological epochs. 




-10 



10 



a, / 1CT 



-500 L 



FIG. 5: The allowed parameter space in this plot is the region 
above the higher solid line and below the lower solid line, with 
h — 0.7 to 2o for a universe containing matter and radiation 
only. The dashed lines show the corresponding result for a 
universe with a non-zero vacuum energy. 



B. Dynamically— coupled theories 



2. Universe containing matter, radiation and a nonzero 
vacuum energy 

A more realistic constraint would involve taking into 
account a late-time period of vacuum domination; so as 
well as p r oc a~ 4 and p m oc a~ 3 we now also have pa — 
constant in the J frame. Hence, 



G eql _ (2uj + 3) 0o _ (2w + 3) 



>eq2 (PO 



G N (2lj + 4) cj) eql (2lu + 4) (j) eql cj) eq2 
_ (2w + 3) 



(2uj + 4) 



(1 + Zegl) (1 + z eq2 ) U+"0(i+*-J 



where we have defined the redshift of matter-radiation 
equality, z eq %, and the redshift of matter-vacuum equal- 
ity, z eq2 , as 



1 + 2, 



eql 



PmO 
PrO 



and 1 + z„ 



„2 



Pad 

PmO 



As above, we still have p r0 = 7.84 x 10 3A g cm 3 but 
now our assumption of spatial flatness gives 



PtotO — PAO + PmO + PrO 



3ffg (4 + 3w)(4 + 2w) 
87tGjv 6(1 + lu) 2 ' 



or, using G N = 6.673 x 10" 11 Nm 2 kg~ 2 and H 
lOO/i kms^Mpc' 1 , 



1 + Zeql = 



f2.39 x 10 4 /t 2 (4 + 3 "l (4 + 2 - ,) -1 



X _|_ PAO 



Taking the value p\o/PmO — 2-7, consistent with WMAP 
observations we the results shown in|SJ 



Constraints on ui(<p) at matter-radiation equality 



Recalling that A 2 
to re-write l|3*3*|) as 



1/0 and a 2 = 2u> + 3 allows us 



the ratio G eq /GN 
and a as 



G = A 2 {l + a 2 ), 

can then be expressed in terms of A 



G e q 

Gn 



A 



eq 



Assuming that In Aq and ln(l + ota) are negligible com- 
pared to lnA eq (i.e. the Universe is close to GR today, 



J) we can write 

In I ^ eq 

\Gn 



2vr 1 



{2uo eq + 3) 



This allows us to constrain ui eq in terms of f3 and a 2 , 
as shown in figure Constraints imposed upon lu and 
a 2 have the same consequences for the evolution of G 
as previously discussed. The parameter (3 controls the 
evolution of w, as can be seen from 1)21 [I. In the limit 
(3 — > it can be seen that u> becomes constant and this 
class of ST theories becomes indistinguishable from the 
BD theory. Constraints upon therefore correspond to 
constraints on the allowed variation of ui and hence G. 

The parameter j3 is taken to be small here so that 
the "kick" on the scalar field during electron-positron 
annihilation can be neglected. These effects have been 
explored by Damour and Pichon in |l8j| , for the case a 2 = 
0, and are expected to have the same result in this more 
general scenario; we will not repeat their analysis of this 
effect here. 
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FIG. 6: The allowed parameter space in this graph is the re- 
gion above the line in the region ui > and below the line 
in the region ui < 0, to 2a. The solid lines correspond to 
a2 — 1CF 11 the dashed lines to 02 = 10~ 13 and the dotted lines 
to 02 = — 10 -11 . Lu eq is the value of u at matter-radiation 
equality and (3 is a parameter of the theory, defined in \21)). 



Wo 








FIG. 7: The allowed parameter space in this graph is the re- 
gion above the line in the region to > and below the line 
in the region u> < 0, to 2a for a universe containing matter 
and radiation only. The solid lines correspond to 02 = 10 -11 , 
the dashed lines to 02 = — 10 -11 ' 523 , and the dot-dashed lines 
to 02 — — 10~ 11155 . The dotted lines show the corresponding 
results for a universe with a non-zero vacuum energy. 



2. Constraints on ujo for a universe containing matter and 
radiation 

The scalar field can be evolved forward in time from 
the time of matter-radiation equality to the present, us- 
ing equation Ij25[l. The two arbitrary constants in this 
expression can be fixed using our limiting values of G eq 
derived above and by assuming that the evolution of 
the scalar field has effectively ceased by this time, i.e. 
%jj' eq = 0, as is the case for the models considered here. 

In order to gain a quantitative limit on vb\ — yJoo, and 
hence on lj q , it is necessary to calculate N eq . From the 
definition of N , we can write 

N eq = - ln(l + z eq ) - \nA eq + In A 

~ - lll(l + Zeq) - ^/3(4>eq ~ V'oo) 2 

where z eq is the red-shift at t eq and we have assumed, as 
before, the term InAo to be negligible. 

It now remains to determine 1 + z eq = p m o / p r o for the 
case co = oj{4>). If we now assume ojq to be moderately 
large, and recall our assumption of spatial flatness, we 
can write 

PtotO — PmO + PrO — ^ — ■ 

For p r0 — 7.84 x I0~ 3i g cm~ 3 this gives p m0 ~ 1.87 x 
1(T 29 h 2 gcm~ 3 . We have 1 + z eq ~ 2.4 x 10 4 h 2 and so 
finally obtain, for h = 0.7, 

Now, evolving ip — vj^ from N eq to No = and using 
we obtain the results shown in figure 



3. Constraints on loq for a universe containing matter, 
radiation and nonzero vacuum energy 

We can repeat the previous analysis for the more real- 
istic case of a universe with a period of late-time vacuum 
domination. We now need the time of matter-radiation 
equality, N eq \ , and the time of matter- vacuum equality, 
N eq2 . 

Assuming spatial flatness, a large u>o, h = 0.7, 
PAo/PmO = 2.7, and p r o as above gives 

N eql ~ -8.06- ^(3{ip eq i -i/>oc) 2 , 

and 

N eq2 ~ -0.33 - i/3(VW2 - ^oo) 2 . 

We can now evolve ip(N) — vb^ from N eql to N eq2 us- 
ing (|25|l and the same boundary conditions as before. 
(In finding N eq2 we used an iterative method to evaluate 
ipeq2 — ipoo)' We then use (|26|l to evolve the field from 
N eq2 to iVo; the constants C and D in l|2t)|l are set by 
matching ip(N) — vj^ and its first derivative with l|25|) 
at N eq2 . Finally, we can calculate luq using (|34|) . The 
results of this procedure are shown in figure 

VI. DISCUSSION 

Using the framework provided by ST theories of grav- 
ity we have investigated the effect of a time-varying G 
during primordial nucleosynthesis. We determined the 
effect on primordial nucleosynthesis numerically, using 
a modified version of the Kawano code [2^, and con- 
strained the parameters of the underlying theory using 
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these results. Our results are consistent with the inter- 
pretation that the abundance of 4 He is primarily only 
sensitive to the value of G at the time when weak inter- 
actions freeze out, for the class of ST models studied. 

Using our numerically determined constraints on the 
evolution of G, we imposed the 2a upper and lower 
bounds shown in figure on the BD parameter to. For a 
constant G (i.e. 02 = 0) we get the bounds w > 332 or 
w < —37, for a universe containing matter and radiation 
only, and the bounds w > 277 or uj < —31, for a uni- 
verse containing matter, radiation and a nonzero vacuum 
energy. As the parameter ai is increased, the strength 
of the upper bound decreases whilst that of the lower 
bound increases; the opposite behaviour occurs as 02 is 
decreased. This is because the strength of the bound on 
uj is essentially due to the allowed value of G at matter- 
radiation equality. If primordial nucleosynthesis allows 
this value to be vastly different from Gn then a signifi- 
cant evolution of G during matter domination, and hence 
a low |w|, is permitted. A value of G eq close to Gn means 
that only a very slow variation of G, and hence large |w| 
is permitted. As 0,2 > corresponds to an increasing G 
during primordial nucleosynthesis this corresponds to a 
higher value of G eq and hence a tighter upper bound on w 
and a looser lower bound (as G decreases during matter 
domination for uj > —3/2 and increases for uj < —3/2); 
a2 < corresponds to G decreasing during radiation 
domination, and so has the opposite effect. 

The more stringent effect of a nonzero value of a-i 
on the upper bounds is due to the current observa- 
tional determinations of the 4 He abundance disfa- 
voring G > Gn during primordial nucleosynthesis (see 
e.g. El). 

The interpretation of the constraints on the more gen- 
eral ST theories is a little more complicated, due to the 
increased complexity of the theories. The bounds on uj(4>) 
at matter-radiation equality, shown in figure are seen 
to be stronger for smaller (3 and weaker for larger for (3 
(assuming f3 to be small enough to safely ignore the ef- 



fect of the e~e + kick analysed by 0])- This should be 
expected asw~ (3~ 2 in the models we are studying. We 
see from the constraints on w at the present day, shown 
in figure [7| that the bounds on w become tighter as j3 
gets very small and large, with an apparent minimum in 
the bound at (3 ~ 0.2. The tight bounds for very small (3 
are due to the tight bounds on uj eq for small (3. The tight 
bounds at large (3 are due to the attraction towards GR 
at late times that occurs for this class of ST theories (see 
e.g. H5). This attractor mechanism is more efficient for 
larger values of (3, as can be seen from pUjl. and at late 
times so w is drawn to a larger value for a larger (3. 

The inclusion of a late-time vacuum-dominated stage 
of the Universe's evolution is to weaken slightly the 
bounds that can be placed on uj at the present day, as 
can be seen from figures |S] and [7\ This weakening of the 
bounds is due to a shortening of the matter-dominated 
period of the Universe's history which is essentially the 
only period, after the effects of the free scalar-dominated 
phase become negligible, during which G evolves. 

We find that the constraints that can be imposed upon 
the present day value of uj from primordial nucleosynthe- 
sis are, for most of the allowed parameter space, con- 
siderably weaker than those obtained from observations 
within the solar system. To date, the tightest constraint 
upon ujq are imposed by Bertotti, less and Tortora 
who find I wq I > 40000, to 2a. This constraint is obtained 
from observations of the Shapiro delay of radio signals 
from the Cassini spacecraft as it passes behind the Sun. 
We consider the constraints imposed upon uj here to be 
complementary to these results as they probe different 
length and time scales, as well as different epochs of the 
Universe's history. 



ACKNOWLEDGEMENTS 

T.C. is supported by PPARC. 



[1] J. K. Webb et al, Phys. Rev Lett. 82, 884 (1999); M. T. 
Murphy et al, MNRAS 327, 1208 (2001); J. K. Webb et 
al, Phys. Rev. Lett. 87 91301 (2001); M. T. Murphy, J. 
K. Webb and V. V. Flambaum, MNRAS 345, 609 (2003). 

[2] It has been suggested (see e.g. |3j) that these observations 
could be due to a different isotopic magnesium abun- 
dance. However, it has also been shown that matching 
the isotopic abundances so as to emulate the apparently 
observed variation in alpha also creates an over abun- 
dance of nitrogen Q . 

[3] T. Ashenfelter, G. J. Mathews and K. A. Olive, Phys. 
Rev. Lett. 92, 41102 (2004). 

[4] Y. Fenner, B. Gibson and M Murphy, MNRAS 358, 468 
(2005). 

[5] P. Jordan, Nature 164, 637 (1949). 

[6] C. Brans and R. H. Dicke, Phys. Rev. 124, 925 (1961). 

[7] P. G. Bergmann, Int. J. Theo. Phys. 1, 25 (1968). 



[8] K. Nordtvedt, Astrophys. J. 161, 1059 (1970). 
[9] R. V. Wagoner, Phys. Rev. D 1, 3209 (1970). 
[10] J. D. Barrow and J. P. Mimoso, Phys. Rev. D 50, 3746 
(1994). 

[11] J. D. Barrow and P. Parsons, Phys. Rev. D 55, 1906 
(1997). 

[12] J. D. Barrow, Phys. Rev. D 47, 5329 (1992). 

[13] J. A. Casas, J. Garcia-Bellido and M. Quiros, Mod. Phys. 
Lett. A 7, 447 (1992). 

[14] J. A. Casas, J. Garcia-Bellido and M. Quiros, Phys. Lett. 
B 278, 94 (1992). 

[15] A. Serna, R. Dominguez-Tenreiro and G. Yepes, Astro- 
phys. J. 391, 433 (1992). 

[16] A. Serna and J. M. Alimi, Phys. Rev. D. 53, 3087 (1996). 

[17] D. I. Santiago, D. Kalligas and R. V. Wagoner, Phys. 
Rev. D 56, 7627 (1997). 

[18] T. Damour and B. Pichon, Phys. Rev. D 59, 123502 



(1999). 

[19] X. Chen, R. J. Scherrer and Gary Steigman, Phys. Rev. 

D 63, 123504 (2001). 
[20] V. Pettorino, C. Baccigalupi and G. Mangano, J. Cosmol. 

Astropart. Phys. 01 014 (2004). 
[21] L. E. Gurevich, A. M. Finkelstein and V. A. Ruban, As- 

trophys. and Sp. Sci. 22, 231 (1973). 
[22] B. Bertotti, L. less and P. Tortora, Nature 425, 374 

(2003) . 

[23] J. D. Barrow, D. Kimberly and J. Magueijo, Class. 

Quant. Grav. 21, 4289 (2004). 
[24] H. Nariai, Prog. Theo. Phys. 40, 49 (1968). 
[25] T. Damour and K. Nordtvedt, Phys. Rev. D 48, 3436 

(1993). 

[26] C. Bambi, M. Giannotti and F. L. Villante, 
|astro-ph/ 0503502 [Phys. Rev. D (to be published)]. 

[27] L. Kawano, Fermilab-pub-92/04-A, 1992. 

[28] D. Kirkman, D. Tytler, N. Suzuki, J. M. O'Meara and 
D. Lubin, Astrophys. J. Suppl. 149, 1 (2003). 

[29] V. Barger, J. P. Kneller, H. S. Lee, D. Marfatia and G. 
Steigman, Phys. Lett. B 566, 8 (2003). 

[30] R. Trotta and S. H. Hansen, Phys. Rev. D 69, 023509 

(2004) . 



11 

[31] C. M. Will, 1993, Theory and Experiment in Gravita- 
tional Physics, 2nd edn. (Cambridge University Press, 
Cambridge) . 

[32] T. Clifton, D. F. Mota and J. D. Barrow, MNRAS 358, 
601 (2005). 

[33] A. R. Liddle, A. Mazumdar and J. D. Barrow, Phys. Rev. 

D 58, 027302 (1998). 
[34] D. J. Fixen et al., Astrophys. J. 473, 576 (1996). 
[35] X. Chen and M. Kamionkowski, Phys. Rev. D 60, 104036 

(1999). 

[36] R. Nagata, T. Chiba and N. Sugiyama, Phys. Rev. D 66, 
103510 (2002). 

[37] O. Zahn and M. Zaldarriaga, Phys. Rev. D 67, 063002 

(2003) . 

[38] R. Nagata, T. Chiba and N. Sugiyama, Phys. Rev. D 69, 
083512 (2004). 

[39] V. Acquaviva, C. Baccigalupi, S. M. Leach, A. R. Liddle 
and F. Perrotta, Phys. Rev. D 71, 104025 (2004). 

[40] C. L. Bennett et al., [WMAP collaboration], Astrophys. 
J. Supplement Series, 148, 1 (2003). 

[41] J. D. Barrow and R. J. Scherrer, Phys. Rev. D 70, 103515 

(2004) . 



